Journal of Proteome Research

Alternative splicing is an important regulatory layer in gene expression, but knowledge on the isoform protein molecules continue to lag their canonical counterparts. An open question is whether alternative protein isoforms feature different half-life than the canonical counterpart, which could indicate differential usage and functional diversification. Here we combined a proteogenomics approach with heavy water-based protein turnover analysis to survey 24 pairs of canonical-alternative protein isoforms in the mouse heart. The results provide a reference on their numerical half-life and also reveal widespread differences in isoform stability.

Ancient DNA (aDNA) has transformed the study of hominin relationships, but its preservation in ancient fossils is often limited. Enamel palaeoproteomics offers an alternative molecular approach for taxonomic analysis. In this study, we re-analyse published DDA mass spectrometry data1 from the Denisovan-attributed Penghu 1 mandible (PXD054412)2 and a Neandertal enamel specimen from Gruta de Oliveira, Portugal (PXD038154) 3. Five AMBN peptides carrying the Valine-273 substitution (V273) were validated in the Penghu 1 enamel, three of which were independently detected in both DDA acquisitions. No V273-containing peptide signal was detected in the Neandertal dataset. Conversely, the ancestral Methionine-273 peptide REDPM[+16]AYG was detected exclusively in the Neandertal specimen. Extracted ion chromatograms, isotopic envelope confirmation, and MS2 fragmentation spectra, all support the reported peptide assignments. AMELY-specific peptides additionally support male sex assignment for both ancient individuals. Together, these results confirm the taxon-specific mutual exclusivity of AMBN V273 and M273 variants across Denisovan and Neandertal lineages, establishing AMBN M273V as a molecularly validated diagnostic marker for Denisovan identification from dental enamel. More broadly, targeted MS1 reanalysis of public proteomics datasets provides a scalable complement to ancient genomics for resolving hominin lineage identity and sex determination when DNA is not preserved

Detecting low-abundance bacterial peptides in samples dominated by human proteins using liquid chromatography tandem mass spectrometry (LC MS/MS) is challenging. Bacterial peptides are typically far less abundant than host-derived peptides, and sequence similarity between bacterial and human peptides can complicate confident attribution to a bacterial origin. To address this problem, we developed and benchmarked an in silico spectral library protein identification workflow for detecting Staphylococcus aureus USA300 peptides against a background of A375 cells. Peptide and protein identifications obtained at a 1% peptide-level false discovery rate were compared between the in silico spectral library search and a conventional sequence database search across a defined bacterial dilution series containing 0, 1%, 2%, 10%, 25%, 50%, and 100% bacterial protein. A human sequence similarity filter was applied to identifications from both search strategies to reduce the likelihood of retaining human derived or human like peptide sequences. Both search strategies showed concentration dependent bacterial detection and shared a practical detection limit of approximately 1 to 2% bacterial protein but exhibited distinct performance profiles. The in silico spectral library search produced fewer bacterial-like background identifications in the human-only control, whereas MS2Rescore substantially increased spectral library identifications at higher bacterial concentrations, exceeding database search at 100% bacterial protein. The partially overlapping peptide sets identified by the two strategies indicate that they provide complementary coverage of the bacterial proteome.

Data-independent acquisition (DIA) mass spectrometry has emerged as a powerful tool for neuropeptidomics, but its success relies heavily on the quality of spectral libraries used for peptide identification. There are inherent challenges to mass spectrometry analysis of crustacean neuropeptides, including the endogenous nature in which they are analyzed, extensive post-translational modification (PTM), and atypical fragmentation patterns. Thus, general-purpose proteomic spectral prediction tools may not perform optimally in the endogenous peptide domain. In this study, we benchmark four widely used spectral prediction platforms, Prosit, MS2PIP, AlphaPeptDeep, and UniSpec, to evaluate their performance in predicting the fragmentation of neuropeptides. Using an empirically derived spectral library from crustacean tissues as reference, we assess model compatibility, dot-product similarity, Pearson correlation, and DIA-based identifications across brain, sinus gland, and pericardial organ samples. Our results reveal that no single model comprehensively captures neuropeptide fragmentation characteristics. While UniSpec showed unexpected strengths due to its inclusion of neutral loss ions, AlphaPeptDeep demonstrated the highest spectral similarity, and MS2PIP and Prosit outperformed in DIA-NN identifications. We further highlight the critical impact of neutral loss fragments, present in over 50% of empirical spectra, and emphasize the need for hybrid spectral libraries that integrate complementary strengths across models. This work provides a foundational framework for optimizing spectral library selection in neuropeptidomics and underscores the importance of model-specific biases when analyzing structurally diverse endogenous peptides.

ShortThe SwissProt database contains a stable 20,418 human protein-coding genes and 42,541 human protein sequences. Ribo-Seq suggests about 7,000 additional, non-canonical Open Reading Frames (ORFs) are present in humans, though only a few of them are confirmed by Mass Spectrometry (MS). Detecting these proteins requires extensive database searches, increasing computational load and inflating False Discovery Rates (FDR). Using the ionbot search engine with the OpenProt database allows for reliable detection of non-canonical proteins while controlling FDR. Ionbot surpasses the Trans-Proteomics Pipeline (TPP) in reproducibility, identifying more peptides and proteins supported by multiple spectra. In addition, open modification searches yield better PSMs compared to closed searches. This work highlights the importance of employing cutting-edge search engines in non-canonical protein research, as well as the value of open modification search in correcting errors in non-canonical protein detection. LongO_ST_ABSBackgroundC_ST_ABSThe SwissProt database reports a quite stable 20,418 human protein-coding genes and 42,541 human protein sequences, figures that have remained stable. New techniques like Ribo-Seq indicate that approximately 7,000 additional, non-canonical Open Reading Frames (ORFs) are translated in humans, few of which have been confirmed by Mass Spectrometry (MS). Detecting these non-canonical proteins requires comprehensive database searches, which increase computational load and False Discovery Rate (FDR). Here, we use the open search engine ionbot in combination with the OpenProt proteogenomics database to reproducibly detect non-canonical proteins while maintaining a well-controlled FDR. ResultsCompared to the current gold standard, the Trans-Proteomics Pipeline (TPP), ionbot shows higher reproducibility, with a higher number of peptides and proteins supported by multiple spectra, and across multiple samples. We observe that PSMs from the open modification search against OpenProt have higher fragment ion intensity correlation compared to PSMs obtained from the closed search, or by only searching canonical proteins. ConclusionsIn this work, we show the potential for open modification searching to correct potential mistakes in non-canonical proteins detection by preventing modified canonical peptides or variants from being incorrectly identified as non-canonical peptides. We also highlight the importance of assessing the FDR of non-canonical identifications separately from canonical ones, as global FDR calculations are biased by the scarcity of non-canonical identifications in each dataset.

The processes such as DNA replication, transcription, and repair are often modulated by specific nuclear proteins, protein-protein interactions (PPIs), and post-translational modifications (PTMs). In Trypanosoma cruzi, however, the nuclear proteome and interactome have not been systematically mapped, limiting the interpretation of nuclear regulatory processes. Here, we report a nuclear proteome resource generated from intact nuclei isolated from T. cruzi and analyzed by high-resolution Orbitrap LC-MS/MS, integrating proteome profiling, computational interaction network inference, and exploratory crosslinking mass spectrometry (XL-MS). Proteome profiling identified 1,734 proteins in the nuclear fraction, including 316 proteins identified with PTM-containing peptides. Subcellular localization prediction and Gene Ontology analysis support nuclear enrichment and highlight functions related to transcription, RNA metabolism, and genome maintenance. The in silico interaction network derived from STRINGDB organizes the proteins into functional clusters, including a histone-associated interaction neighborhood. In parallel, XL-MS identified 26 residue-resolved interprotein crosslinks involving 36 proteins and detected PTMs at or near linked residues. Together, these data support reuse for comparative nuclear proteomics, multi-omics integration, and prioritization of candidates for future functional studies.

Affinity purification-mass spectrometry (AP-MS) enables the characterization of protein-protein interactions (PPIs), and the ease and sensitivity of such experiments has progressively increased. Beyond steady-state interactions of target proteins, a strong interest has emerged in monitoring how PPIs change upon significant biological perturbations, such as in disease contexts or small molecule modulation of the target protein. These perturbations likely not only induce PPI changes but can also lead to altered expression of proteins not of direct interest. Changes in protein abundance may alter which proteins adsorb to the affinity purification matrix, and due to the sensitivity of modern mass spectrometers, these differential ''background binders'' can masquerade as differential interactors. Contemporary approaches often do not account for differences in the background proteome, potentially inflating the number of false positives and negatives reported. Here, we provide technical considerations for the reliable annotation of dynamic PPIs, using the O-GlcNAc transferase (OGT) as a case study. We describe the installation of affinity epitope tags on endogenous OGT in mouse embryonic stem cells (mESCs), which we then apply for OGT interactor identification via AP-MS. We show that accurate representation of the bead background, which depends on the affinity matrix in use, is critical for elimination of false positive and false negative PPIs. This became even more pertinent as OGT PPI dynamics were measured under OGT catalytic inhibition via OSMI-4, which is known to perturb gene expression. The proteomes of OSMI-4-treated and control-treated mESCs differed, leading to distinct bead backgrounds in which the differential background proteins appeared as interaction gains or losses. These false positives were resolved by incorporating straightforward experimental controls through a practical statistical framework, allowing for a direct and confident comparison between treatment conditions. Incorporating these considerations into workflows investigating PPI dynamics will improve data fidelity and reproducibility.

As the most prevalent neurodegenerative disorder worldwide, Alzheimer's disease (AD) remains incompletely understood at the proteome level. In particular, current studies on amino acid (AA) substitutions have predominantly relied on genomic and transcriptomic profiling. Proteome-scale AA substitutions remain largely uncharacterized. Nevertheless, alterations at the DNA and RNA levels cannot fully recapitulate the spectrum of AA substitutions observed at the proteome level. This critical research gap persists largely due to the inherent analytical challenges posed by large-scale proteomic datasets. In this study, we address this limitation by analyzing two independent AD proteomic datasets, AMP-AD and PXD013753, using PIPI-C, an open-search mass spectrometry engine capable of resolving multiple co-occurring modifications per peptide. We introduce a pipeline that enables the characterization of AA substitutions at the proteome level and the dissection of regulatory functions of key proteins with such substitutions. In both datasets, we observe that, after controlling for ambiguous post-translational modification mass shifts, the N>M substitution is the most frequent variant among N>X substitutions in the AD data. Furthermore, among the 701 overlapping proteins shared by the two datasets, we identified 6 literature-reported substitution sites, including residues 242 and 352 in glial fibrillary acidic protein, site 370 in actin gamma 1, as well as sites 111, 115, and 116 in hemoglobin subunit beta. Since most reported AA substitutions rely on genomic or transcriptomic evidence, and our pipeline adopts rigorous false-positive controls, those unreported substitutions are presumably detectable only via proteomic methods, emphasizing the unique value of our proteome-based substitutomics workflow for AD research.

Proteolysis can be completed in a single-step, rapid (5-30 minutes) reaction using hyperthermoacidic archaeal (HTA) proteases. Protocol simplicity and stability of HTA proteases in ambient storage conditions removes many logistical barriers to widespread deployment of proteomic analysis. HTA protease cleavage selectivity is complicated but quantitatively reproducible and generates peptides that negligibly overlap with trypsin. Tryptic and HTA protease peptides emphasize different proteins and provide complementary coverage of amino acid sequence and post-translational modification (PTM) sites. Digestion of a sample's separate aliquot with HTA protease is a minor addition to trypsin-based workflows, but LC-MS analysis of both sets of peptides is not practical for large cohorts. This manuscript evaluates the identifications and reproducibility of LC-MS analysis of tryptic and HTA protease peptide mixtures. Adding HTA protease peptides to 20% increased total protein identifications by 24% on average in neat plasma and identified 216 additional non-artifact modifications compared to conventional trypsin-only bottom-up proteomic analysis. This approach is reproducible and does not increase sample preparation or analysis time. Peptides from each protease present independent opportunities to quantify overlapping proteins and a quantification model was developed in which the protease with the most complete data is selected for each protein while the second protease is used to infer missing values. This cross-protease inference is more accurate than the alternative minimum, mean, or nearest neighbor imputations. Finally, biomarkers discovered through HTA protease peptides are amenable to rapid HTA protease only proteomic workflows thus streamlining translation from discovery to fieldable applications.

SummaryContemporary proteomics methods can now generate large-scale DIA datasets of thousands of files that demand substantial computational resources for efficient analysis. Spectronaut is a widely used platform for DIA data processing; however, large-scale searches are often constrained by computational performance and long execution times when run on single workstations. Here, we present Spectronaut-nf, a Nextflow-based pipeline that enables scalable and parallelized execution of Spectronaut analyses across high-performance computing (HPC) environments. The workflow divides directDIA analysis into modular stages, including spectral library generation, DIA searching, and merging results, allowing efficient distribution of tasks across multiple compute nodes. Benchmarking using 72 diaPASEF raw files using typical hardware demonstrated that Spectronaut-nf completed searches in 23.77 hours, compared with 39.09 hours on a Windows workstation and 67.04 hours on a single-node Linux HPC setup. Stress testing with 1,037 diaPASEF raw files further demonstrated the scalability and robustness of the workflow for large proteomics datasets. Across platforms, protein and peptide identifications remained consistent, with only minimal variability attributable to platform-specific differences. Overall, Spectronaut-nf provides a flexible, scalable, and efficient framework for high-throughput DIA proteomics analysis in HPC environments. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=129 SRC="FIGDIR/small/741433v1_ufig1.gif" ALT="Figure 1"> View larger version (33K): org.highwire.dtl.DTLVardef@4a7107org.highwire.dtl.DTLVardef@14287d8org.highwire.dtl.DTLVardef@e484e9org.highwire.dtl.DTLVardef@d21c07_HPS_FORMAT_FIGEXP M_FIG C_FIG

Quantitative proteomics requires decisions about normalization, missing-data handling, and statistical modeling that can change the reported results. ProteoScopeR is an R package and Shiny application that connects these decisions in a traceable workflow and complements downstream exploration in xOmicsShiny. Side-by-side comparisons expose changes in distributions, feature retention, effect estimates, and selected protein sets. An optional external artificial-intelligence assistant reviews exported evidence and proposes settings for researcher approval; statistical calculations remain in R. A group-only case study used DIA-NN-derived aqueous humor data from 69 samples and 3,667 imported protein groups. Six normalization methods and seven missing-data strategies were screened, and the approved cyclic-loess analysis compared five differential workflows. For nAMD versus pmCNV, limma and proDA selected 54 and 37 protein groups, respectively, with 31 shared selections at adjusted P < 0.05 and absolute log2 fold change [&ge;] 0.5. Input routing, scientist decisions, locked settings, and an execution receipt preserve the analysis record. These comparisons describe sensitivity to analytical choices rather than identify a universally superior method.

Mass spectrometry (MS) is the leading technology for identifying proteins in complex biological samples. It relies on the use of tandem MS alongside a reference database of canonical protein sequences to computationally identify peptides and their parent proteins. The canonical sequences represent the most widely expressed and functionally validated forms of proteins. Consequently, disease-induced or disease-supportive variants, such as those associated with cancer, will evade detection if they are absent from the database. To address this challenge, this study introduces a revised release of the Unkown Mutation Analysis (XMAn) database by incorporating coding missense and nonsense mutations from the latest versions (v103) of the COSMIC Genome Screen Mutants (GSM) and Cancer Gene Census (CGC) datasets in two distinct FASTA-formatted peptide databases comprising 3,848,499 and 312,658 variants, respectively. The mutated peptides were matched to reviewed, non-redundant UniProt Homo sapiens protein entries (18,362 and 746), and characterized in terms of nucleotide- and amino acid mutation frequencies, peptide length distributions, and associations between specific single-nucleotide (SNV) and single amino acid (SAAVs) variants. Applied to the analysis of MDA-MB-231 breast cancer cell-membrane protein fractions, the database enabled the identification of 300+ high-quality variant peptides - several localized to functional protein-binding and catalytic domains - and 23 aberrant protein products mapped to the CGC dataset. The database is hosted and available for download on Zenodo (XMAn/gsm doi: 10.5281/zenodo.21781023; XMAn/cgc doi: 10.5281/zenodo.21781514) or can be accessed through https://sites.google.com/vt.edu/xman-db/home.

Plasma proteomics is increasingly used for biomarker discovery and predictive modeling, yet diurnal protein trajectories remain insufficiently characterized. In our review of recent proteomic biomarker studies, 43% of the identified biomarkers had previously been reported to display 24-h rhythmicity. We demonstrate that ignoring these short-term dynamic effects compromises the robustness of reported models predicting health outcomes. We integrated a population-scale multi-ethnic longitudinal cohort with repeated measures over 10 years, with two cohorts of healthy adults undergoing frequent plasma sampling across days under controlled circadian, sleep and food-intake conditions. This design enabled estimation of short-term intraindividual variability (ST), long-term intraindividual variability (LT), population-level variability (POP) and genetic effects (GEN) across 7,289 protein targets. ST, LT, POP, and GEN define diverse protein trajectories, including rapid dynamics, long-term change, and individual-specific signatures. Using and generalizing this framework will facilitate covariate selection, study design, biomarker prioritization, and variability-aware modeling by users of proteomic data. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=92 SRC="FIGDIR/small/741635v1_ufig1.gif" ALT="Figure 1"> View larger version (26K): org.highwire.dtl.DTLVardef@55520eorg.highwire.dtl.DTLVardef@17e4bacorg.highwire.dtl.DTLVardef@9a0d0borg.highwire.dtl.DTLVardef@1ce88fb_HPS_FORMAT_FIGEXP M_FIG C_FIG

Proteomics studies increasingly use different technologies, including mass spectrometry, the antibody-based Olink platform, and the aptamer-based SomaScan platform. However, downstream analysis often depends on platform-specific scripts or point-and-click tools, limits flexibility and makes analyses difficult to reuse across platforms. We developed ProtPipe2, a downstream proteomics analysis framework for mass spectrometry, Olink, and SomaScan data that is available as both an R package and an interactive web application. ProtPipe2 supports quality control, preprocessing, imputation, batch correction, statistical analysis, visualization, pathway analysis, and result export. ProtPipe2 expands the original Protpipe workflow to additional proteomics platforms and provides reusable R functions for imputation, batch correction, and other downstream analyses. Its redesigned web application provides user-controlled methods and parameters through the same R-based workflow. Case studies across mass spectrometry, Olink, and SomaScan datasets demonstrated a consistent downstream analysis workflow across platforms. By supporting three proteomics platforms through both graphical and scripted workflows, ProtPipe2 makes downstream analysis easier to perform, repeat, and review. ProtPipe2 is freely available at https://nih-card.github.io/ProtPipe2/.

There are many databases of post-translational modifications, but each has different protein and PTM coverage, making the current state of post-translational modifications difficult to obtain, assess, or harmonize. ProteomeScout has been a relational database that harmonized a myriad of resources to address this challenge. In this update, we wished to create a lightweight, sustainable open source resource for the community that retained the power of harmonization, protein visualization, and Python tools for interacting with the database. ProteomeScout-v4 presents an ecosystem of open-source tools for aggregating and updating proteins and post-translational modifications and presenting the web-based and code-based access tools for interacting with the data. ProteomeScout has moved from using a complex relational database to a flat-text file format. The database file update is version controlled on FigShare and all features have processes for continuous updating, including migrating PTMs to updated protein records, when external databases are updated. The open-source nature and flexibility allow for user-based control and community-based feature integration. Fast and easy analysis is possible with the API, with example analyses provided here.

Serum proteomics is increasingly used for minimally invasive biomarker discovery and disease phenotyping, and the choice of serum preprocessing workflow can shape proteome depth, quantitative characteristics, and downstream biological readouts. However, disease-oriented comparisons within a single cohort remain limited. Here, we compared four serum preprocessing workflows--Top14 depletion (TOP14D), tomato lectin affinity purification (TomAP), and two nanoparticle-based enrichment workflows (NPA and NPB)--using serum from six patients with systemic juvenile idiopathic arthritis (sJIA) and six age- and sex-matched healthy controls, and analyzed them using unified data-independent acquisition mass spectrometry (DIA-MS) and a statistical pipeline. We evaluated proteome depth, missingness, quantitative characteristics, group separation, differential abundance signatures, pathway enrichment, curated sJIA-related gene set coverage, pre-ranked gene set enrichment analysis (GSEA) results, and detection of inflammasome/interferon-related proteins. TomAP yielded the greatest proteome depth (7612 proteins), followed by NPB (6735 proteins) and NPA (6602 proteins), whereas TOP14D yielded the smallest protein set (3303 proteins). Principal component analysis (PCA) showed a separation between the sJIA and control groups for all workflows. Differentially expressed proteins (DEPs) showed limited overlap, with only 75 DEPs common to all four workflows. Functional enrichment patterns were workflow-dependent; TOP14D and TomAP mainly captured neutrophil/myeloid and inflammatory processes, whereas NPA and NPB captured RNA processing- and translation-related signals. TomAP showed relatively broad coverage and positive enrichment of curated sJIA-related gene sets associated with inflammation, innate immunity, and macrophage activation syndrome (MAS). Inflammasome/interferon-related proteins, including NLRC4, PYCARD, GSDMD, MEFV, IL-18, OAS3, and MYD88, showed workflow-dependent detectability and differential abundance. These findings support a disease-oriented benchmark for fit-for-purpose workflow selection according to the disease axis and analytical objective rather than proteome depth alone.

Mass spectrometry-based chemical proteomics is a powerful method to analyze proteins labelled by small molecules to identify protein targets of active compounds and to profile protein post-translational modifications. The throughput and high protein input for chemical proteomics workflows has been often a limiting factor for application of the technology for specialized and difficult to culture cell lines. The high protein input was necessary to gain significant difference of noise to signal ratio in proteomics readout. Here, we describe a general chemical proteomics workflow, which is performed in 96-well plate and necessitate only 25 g of protein input to profile post-translationally modified proteins including abundant O-GlcNAcylated proteins as well as low abundant AMPylated proteins. The workflow integrates advances in Cu(I)-catalyzed azide-alkyne cycloaddition to minimize chemical side-reactivity of the click reaction and data-independent acquisition mode during LC-MS/MS measurement. An iterative optimization of protein clean-up on carboxylate-coated paramagnetic beads led to significant saving of the beads usage and lowers the unspecific protein background that resulted in sensitivity gain.

Multicenter studies are essential for benchmarking analytical workflows, yet their interpretation is often confounded by the combined effects of experimental protocols and instrumentation. To address this challenge, we introduce a simple normalization-based analytical framework, the recovery metric ({rho}), designed to decouple protocol driven effects from instrument dependent variability. We applied this framework to the 13th Proteomics Multicentric Experiment (PME13), a large multicentric proteomics dataset generated across 27 laboratories using high sensitivity workflows and varying sample preparation protocols. By leveraging a common digested reference sample, {rho} enables direct cross-comparison of all datasets on a unified scale, effectively minimizing instrument-related biases. Using this approach, we demonstrate that apparent instrument dependent trends are largely removed when evaluated through {rho}, revealing consistent protocol driven effects across laboratories. Statistical modeling identified key variables influencing {rho}, including sample input amount, reduction and alkylation, and the use of n-dodecyl-{beta}-D-maltoside (DDM). While DDM was associated with improved {rho}, reduction and alkylation and additional handling steps led to reduced performance, particularly at low input levels. We further highlight practical considerations for the application of ratio based normalization, including the occurrence of values exceeding theoretical bounds, which reflect deviations from underlying assumptions and require appropriate filtering. Overall, this work establishes a generalizable analytical strategy for disentangling confounding factors in multicentric datasets and provides practical guidelines for optimizing high sensitivity proteomics (HSP) workflows. The proposed framework is broadly applicable to other analytical fields where cross laboratory comparability is required.

Proteomics is a powerful tool for profiling the protein landscape underlying cellular and physiological processes. Advances in mass spectrometry-based methods, such as tandem-mass-tags (TMT) isobaric labeling-based proteomics at both MS2 and MS3, label-free quantification (LFQ) with data-dependent-acquisition (DDA) or data-independent acquisition (DIA)-based proteomics, and stable isotope labeling with amino acid (SILAC) spike-in-assisted proteomics, enable quantitative proteome comparisons across biological conditions. These five quantitative proteomic methods offer distinct advantages in coverage, throughput, and quantitative reproducibility. Here, we treated them as orthogonal approaches, retaining protein candidates with significant, consistent-direction changes supported across multiple independently instrumented methods to identify high-confidence protein-level changes associated with longevity. Insulin/insulin-like growth factor signaling (IIS) is a highly conserved longevity regulatory pathway from Caenorhabditis elegans to humans. By quantitatively profiling IIS mutants in C. elegans and performing genetic screening, we obtained a systems-view of proteomic changes linked to longevity and identified new longevity regulators. This systematic cross-method comparison also serves as a valuable resource for guiding the selection of quantitative proteomics approaches.

Mass spectrometry-based proteomics increasingly supports study designs with repeated measurements, such as longitudinal sampling and spatial tissue profiling. These designs offer rich biological insights but introduce statistical challenges due to within-subject correlations and variance heterogeneity across time points or anatomical locations. Despite the growing use, best practices for analyzing repeated-measures proteomics data remain limited. We systematically compare commonly applied methods, including paired t-tests and linear mixed-effects models with correlation induced by random effects, to more flexible covariance-pattern models with unstructured covariance. Through simulation studies, we evaluate false-positive rates and statistical power under varying conditions of covariance heterogeneity. Our results demonstrate substantial limitations of conventional approaches and highlight the robustness and superior performance of unstructured covariance pattern models for repeated measures proteomics. This approach accommodates diverse study designs supporting accurate inference and reproducibility. We demonstrate that conventional F-tests can be severely biased, not just for small sample sizes but also moderately large ones, while using pairwise t-tests or applying parametric bootstrapping to evaluate the p-value of the F-tests mitigates the problem. We provide real world data examples showing how to perform valid and informative statistical analyses for temporal and spatial proteomics studies with repeated measurements.